Generation and control of Greenberger-Horne-Zeilinger entanglement in superconducting circuits 
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Going beyond the entanglement of microscopic objects (such as photons, spins, and ions), here we propose 
an efficient approach to produce and control the quantum entanglement of three macroscopic coupled super- 
conducting qubits. By conditionally rotating, one by one, selected Josephson charge qubits, we show that their 
Greenberger-Horne-Zeilinger (GHZ) entangled states can be deterministically generated. The existence of GHZ 
correlations between these qubits could be experimentally demonstrated by effective single-qubit operations fol- 
lowed by high-fidelity single-shot readouts. The possibility of using the prepared GHZ correlations to test the 
macroscopic conflict between the noncommutativity of quantum mechanics and the commutativity of classical 
physics is also discussed. 



PACS number(s): 03.65.Ud, 03.67.Lx, 85.35.Ds 

Introduction. — Entanglement is one of the most essential 
features of quantum mechanics and has no analogue in classi- 
cal physics. Mathematically, it means that the wave function 
of a system composed of many particles cannot be separated 
into independent wave functions, one for each particle. Physi- 
cally, entangled particles can display remarkable and counter- 
intuitive quantum effects. For example, a measurement made 
on one particle collapses the total wavefunction and thus in- 
stantaneously determines the states of the other particles, even 
if they are far apart. 

The existence of entanglement has been experimentally 
demonstrated |1] with, e.g., two photons separated far apart 
(e.g., up to 500 m) and two closely-spaced trapped ions (e.g., 
separated a few micrometers apart). The obvious violation of 
Bell's inequality in these two-qubit experiments statistically 
verifies the conflict between the locality of classical physics 
and the non-locality of quantum mechanics. Only recently, 
the experimental study of entanglement has been successfully 
extended to a system composed of more than two qubits. For 
example, three-photon Greenberger-Horne-Zeilinger (GHZ) 
entangled states |2] have been demonstrated, and then used 
to test the conflict between classical local-realism and quan- 
tum non-locality using definite predictions 1 3], rather than the 
statistical ones based on Bell's inequalities. Yet, besides the 
problem of detector efficiency, the expected GHZ state in op- 
tical experiments could not be deterministically prepared |2] 
because: i) each entangled photon-pair was generated in a 
small subset of all pairs created in certain spontaneous pro- 
cesses, and ii) the nondeterministic detection of a trigger pho- 
ton among two pairs of entangled photons was required. 

Instead of fast-escaping photons, massive or macroscopic 
quantum systems |4] have also been extensively studied to 
realize controllable multipartite quantum entanglement. The 
three-qubit entanglement of microscopic Rydberg atoms |5] 
and trapped-ions |6] was prepared experimentally. Moreover, 
the GHZ state of massive macroscopic "particles" has also 
been demonstrated in liquid NMR 1 7 ] . However, the existence 
of nonlocal correlations in these "particles" cannot be settled, 
as the correlated information between them will be completely 



mangled in their readouts of ensemble averages. 

Superconducting qubits |8] provides an attractive platform 
to control the genuine (rather than ensemble-pseudo-pure) 
macroscopic quantum state. The sizes of the present "par- 
ticles", e.g., Cooper-pairs boxes, and the distance between 
them are typically on the order of microns. If the interbit cou- 
plings are switchable, then methods OlSlal- working well in 
photon- and trapped-ion systems, could be applied |9] to gen- 
erate and verify the GHZ entanglement between the Joseph- 
son qubits. However, in all published (so far) experiments the 
interactions between Josephson qubits 1 8] are fixed (either ca- 
pacitively or inductively), and thus the usually required single- 
qubit gates cannot, in principle, be strictly implemented. 

For the currently-existing experimental circuits with 
always-on coupling, here we propose an effective approach 
to deterministically generate three-qubit GHZ states by con- 
ditionally rotating the selected qubits one by one. The exis- 
tence of the desirable GHZ entanglement is then reliably ver- 
ified by using effective single-qubit operations. The prepared 
GHZ entanglement should allow to test quantum nonlocality 
by definite predictions at a macroscopic level. 

Preparation of GHZ states. — We consider the three-qubit 
circuit sketched in Fig. 1 ; that is^ only adding one qubit to 
the experimentally-existing one 1 10]. Three superconducting- 
quantum-interference-device (SQUID) loops with control- 
lable Josephson energies produce three Josephson qubits, fab- 
ricated a small distance apart (e.g., up to a few microme- 
ters lHoll . as the case of entangled trapped ions in Ref. JfJ) 
and coupled via the capacitances Cn and C23. The dynam- 
ics of the system can be effectively restricted to the subspace 
spanned by the computational basis, and be thus described by 
the following simplified Hamiltonian 
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•> 0.5, is the effective charging en- 
ergy of the jth qubit, whose effective Josephson energy is 
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E^p =2e K j > cos(7r<i>j/$o) with £j j the Josephson energy of 
the single-junction and $ the flux quantum. The effective 
coupling energy between the jth qubit and the (j + l)th 
one is Kj,j + i=e 2 Cj^ +l . Above, Cs j is the sum of all ca- 
pacitances connected to the jth box, and other effective ca- 
pacitances are defined by Cs 1 = + Cf 2 Cs 3 /C), 
Cs 2 = C , /(Cs 1 Cs 3 ), Cs 3 = Cs 3 /(l + C| 3 Cs 1 /C), Cyi — 

C/(Ce 3 Ci2), C23 = C/(C-E 1 C 2 3), C13 = C/(Ci2C23), 

with C=n- = i^s J - C\ 2 C^ - C 2 2 3 C7 Sl . The pesu- 
dospin operators are defined as crS J ' ) = | ) (Oj | — and 
(Tx — |0j) (lj I + 1 Ij) (0j\. As the interbit-couplings are always 
on, the charge energy Eq of the jth qubit depends not only 
on the gate-voltage applied to the jth qubit, but also on those 
applied to the other two Cooper-pair boxes. Compared to 
the coupling Kj j +1 between nearest-neighboring qubits, the 
interaction of two non-nearest-neighbor qubits (i.e., K13 — 
e 2 /Ci3 between the first and the third qubits), is very weak 
and thus has been safely neglected [11]. Indeed, for the typical 
experimental parameters: Cj ~ 600 aF, C m ~ 30 aF, and 
C g = 0.6 aF in Ref. Gil, we have K 13 / K 12 =K 13 / K 23 < 
C m /Cj = 0.05 and K 12 /2sj ~ 1/4. 

In principle, the coupled qubits cannot be individually 
manipulated, as the nearest-neighbor capacitive couplings 
Kjj + i are sufficiently strong. However, once the state of the 
circuit is known, it is still possible to design certain operations 
for only evolving the selected qubits and keeping the remain- 
ing ones unchanged. Our preparation begins with the ground 
state of the circuit |^(0)) = |000), which can be easily ini- 
tialized. The expected GHZ state could be produced by the 
following simple three-step pulse process 
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The first evolution U 2 (t 2 ), with sin[£^ 2) t 2 / \2h)]=±l / V2, 
is used to superpose two logic states of the second qubit. 

This is achieved by simply using a pulse that switches on 

(2) 

the Josephson energy Ej and sets the charging energy 
E^ ] =-2(K 12 +K 23 ). The second (or third) evolution Ui{t\) 
(or il 3 (t 3 )) is achieved by switching on the Josephson en- 
ergy of the first (third) qubit and setting its charging en- 
ergy as Eq =2K\ 2 (or Eq — 2K 23 ). The corresponding 
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duration is set to satisfy the conditions wo[Ej tj/{2K)]- 

and cos(7jtj-/a)=l, with 7^- = ^ (2K j2 ) 2 + (E { ] ] /2) 2 , with 
j = 1, 3, in order to conditionally flip the jth qubit; that is, 
flip it if the second qubit is in the |1) state, and keep it un- 
changed if the second qubit is in the |0) state. 

The fidelity of the GHZ state prepared above can be exper- 
imentally measured by quantum-state tomography lla. 171 fl2ll . 
However, it would be desirable to confirm the existence of 
a GHZ state without using tomographic measurements on a 




FIG. 1: (Color online) Three capacitively-coupled SQUID-based 
charge qubits. The quantum states of three Cooper-pair boxes (i.e., 
qubits) are manipulated by controlling the applied gate-voltages 
Vj (j = 1,2,3), and external magnetic-fluxes $j (penetrating the 
SQUID loops). The circuit can be generalized to include more qubits. 



sufficient number of identically prepared copies. Optical ex- 
periments |01 have achieved this via single-shot readout and 
we propose a superconducting-qubit analog of this approach. 
The single-shot readout of a Josephson-charge qubit has been 
experimentally demonstrated 1 1 3 ] by using a single-electron 
transistor (SET) [14]. Before and after the readout, the SET 
is physically decoupled from the qubit. The GHZ state gener- 
ated above implies that the three SETs, if they are individually 
coupled to each one of the three Cooper-pair boxes at the same 
time, will simultaneously either receive charge signals or re- 
ceive no signal. The former case indicates that the circuit is 
in the state | 111), while the latter one corresponds to the state 
1 000). However, the existence of these two terms, 1 1 1 1) and 
1 000), in these single-shot readouts, is just a necessary but not 
yet sufficient condition for demonstrating the GHZ entangle- 
ment. Indeed, a statistical mixture of those two states may also 
give the same measurement results. In order to confirm that 
the state (2), e.g., \ipQ HZ ), is indeed in a coherent superposi- 
tion of the states 1 000) and | 111), we consider the following 
operational sequence 
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which is similar to the verification of the optical GHZ correla- 
tions H. Above, P 2 = |1 2 )(1 2 | is a projective measurement 
of the second qubit. The suffixes are introduced in the sec- 
ond and third steps to denote the order of the qubits. When 
we finally readout the first and third qubits at the same time, 
the simultaneous absence of the terms IO1O3) and 1 1 1 13) due 
to destructive interference indicates the desired coherent su- 
perposition of the terms in the prepared GHZ state (2). The 
question now is how to realize the required single-qubit oper- 
ations Uj = exp[iTTai^ /4], j — 1, 2, 3, keeping the remain- 
ing qubits unchanged, in this circuit with untunable interbit 
interactions (like the currently available experimental ones). 

In order to effectively implement the single-qubit rotation 
U 2 performed only on the second qubit, while keeping the first 
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and third qubits unchanged, we let the circuit evolve under the eigenvalue +1. Thus, classical reality implies that 



Hamiltonian H 2 = -efa^ + K 12 o z 1] a z 2) + K 23 af ] af\ 
by only switching on the Josephson energy of the second 

qubit, e.g., Ef = 2ef. Since C12 = K 12 /{2ef) < 

1, C23 = K 2i /{2e { j ] ) < 1 (e.g., < 1/4 for the typical ex- 
perimental parameters 1 10]), we can treat the second and third 
terms in H 2 as perturbations of the first one there. Indeed, 
neglecting quantities smaller than the second-order perturba- 
tions 1 1511. th e Hamiltonian H 2 can be effectively approxi- 
mated \o\\(& 
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In the state (2) the logic states of the first and third qubits 
are always identical. Thus, by setting the corresponding dura- 
tion r 2 as r 2 = hir/{4e ( ? [l + 2( 2 12 + 2C| 3 + 4(12(23] }, the 
required single-qubit operation U2 — exp(— ziT^^/fr) — 

(2.) 

exp(iTT(Tx /4) could be effectively performed on the sec- 
ond qubit in state (2). Similarly, the Hamiltonian H13 = 



Ej=i a{~ £ j ux + Kj 2 a y i' induced by simultane- 
ously switching on the Josephson energies of the first and third 
qubits, can be effectively approximated to 
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by neglecting the higher-order terms of Q 2 = Kj 2 /{2s*p) < 
1, with j = 1, 3. The shifts of Josephson energies AEj — 

4£j \f 2 cr z 2 ^ depend on the state of the second Cooper- 
pair box, which collapsed into the state |0) after the pro- 
jective measurement P 2 — 1 12} {I2 1 (because such a mea- 
surement tunnels the existing excess Cooper-pairs into the 

" (13) 

connected SET). Thus, the effective Hamiltonian H^ s 
yields the evolution J7i3(ti3J = exp(— iH cS T13/H) — 
n j=1 ^exp{zr 13 [£y ) (l + 2C| 2 )]CTi j) /^}- Obviously, if the 
duration T13 satisfies the condition txs[ej(1 + 2(| 2 )]//i — 
7r/4, then the required single-qubit operations Uj — 
exp[i7T<7x /4] could be simultaneously implemented. 

Possible application. — The prepared GHZ state, e.g., 
\ipQ HZ ), should allow, at least in principle, to test the macro- 
scopic conflict between the noncommutativity of quantum 
mechanics and the commutativity of classical physics by def- 
inite predictions |3|. Using the EPR's reality criterion, each 
observable corresponds to an "element of reality" (even if it is 
not measured). That is, the quantum operators crtP , (a = 
x,y,z;j — 1,2,3) are linked to the classical numbers 

(i) (?) 
rria , which have the value +1 or —1. The so called a a , - 

measurement is the projection of the quantum state into one of 
the eigenstates of a a ■ The prepared GHZ state is the eigen- 



state of the three operators: A y 



(1) (2) (3) , . 

cr x cf v oh , and A xxy 



_.(i)_.(2)_(3) A 

Uy Ux Ux 1 ^xyx 



(1) (2) (3) 

cr x Ox <J V , with a common 



1 = (m y ^m^m^)(mPm^m^)(m^m^m^) 

The second formula indicates that, if we perform the 
Ilj=i a v -measurement (i.e., yyy-experimenf) on the state 
\i>QH Z }, the eigenstate |— ) only shows in pairs. Here, |+) 
(or |— )) denotes the eigenstate of the operator a y with eigen- 
value +1 (or —1) and corresponds to the classical number 
m y = +1 (or —1). While, for this yyy-experiment quantum- 
mechanics predicts that the state |— ) never shows simultane- 
ously in pairs, because the prepared GHZ state can be rewrit- 
ten as \^+ HZ ) = (|++-) + |+-+) + |-++) + |---»/2. 
Obviously, this contradiction comes from the fact that the ob- 
servable o~x^ anti-commutes with the observable Oy and the 
operator identity 

which is "opposite" to its classical counterpart. 

The protocol described above could be directly (e.g., for 
the optical system yl) performed by reading out the eigen- 
states of the operators a x and a y , respectively. However, in 
the present solid-state qubit, the eigenstates of a z are usu- 
ally read out. Thus, additional operations, e.g., the Hadamard 
transformation S x = (o z + a x ) / v2, and the unitary transfor- 
mation S y = [(1 + i)I + (1 — i) 53 fa]/ (2V2), are required 
to transform the eigenstates of a x and a y to those of cr z , re- 
spectively. These additional single-qubit operations could be 
implemented by combining the rotations of the selected qubit 
along the x-axis (by using the effective Hamiltonian proposed 
above) and those along the z-axis (by effectively refocusing 
the fixed-interactions 1 15]). 

Conclusion and Discussions. — The experimental realiza- 
tion of our proposal for producing and testing GHZ correla- 
tions is possible, although it may also face various technolog- 
ical challenges, like other theoretical designs [17] for quantum 
engineering. Of course, the fabrication of the proposed circuit 
is not difficult, as it only adds one qubit to experimentally- 
existing superconducting nanocircuits [10]. Moreover, rapidly 
switching on/off the Josephson energy, to realize the fast 
quantum manipulations, is experimentally possible. In fact, 
assuming a SQUID loop size of 10 (^m) 2 , changing the flux 
by about a half of a flux quantum in 10~ 10 s, requires sweeping 
the magnetic field at a rate of 10 5 Tesla/s, almost reachable by 
current techniques II 811. 

Also, the prepared GHZ states are the eigenstates of the 
idle circuit (i.e., no operations on it) without any charge- 
and Josephson energies (by setting the controllable param- 
eters as §j = $o/2 and n g . = 0.5) and thus are rela- 
tively long-lived, at least theoretically. Indeed, the couplings 

V = Y^=i a ^ (Y^—i ^jkXkj between the relevant baths 
and the circuit commute with the non-fluctuating Hamiltonian 
of the idle circuit H = J2j=i 2 Kj.j+i o z ^ a{^ +1 \ Here, 
Xjk equals to either 1 for j — k or C^ k Cjk for j ^ k, and 
x k = {eC gk /C Sk )J2u Jk (9u k d 'i k + gu k a^ k ) with a Wk ,al k 
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being the Boson operators of the kth bath, and g UJk the cou- 
pling strength between the oscillator of frequency u>^ and 
the non-dissipative system. Thus, only pure dephasing, i.e., 
the zero frequency value of the noise spectrum contributes 
to overall decoherence rates I19B . However, the working fre- 
quency of the present circuit is always non-zero. This implies 
that the lifetimes of the prepared GHZ correlations are still 
sufficiently long, and thus various required quantum manipu- 
lations could still be coherently implemented. 

Perhaps, the biggest challenge comes from the fast single- 
shot readouts 113112(1 of multi-qubits at the same time. This is 
a common required task of almost all quantum algorithms and 
an important goal for almost all physical realizations of quan- 
tum computing. In order to avoid the crosstalk between qubits 
during the readouts, the readout time t m should be "much" 
shorter than the characteristic time t c ~ H/Kj^+i of commu- 
nications. This requirement has been achieved by the existing 
phase-qubit circuits EUl : t m ~ Ins, and t c ~ 4ns for the 
demonstrated coupling energy K ~ 80 MHz. For the existing 



charge-qubit circuits 1 10], where the interbit coupling-energy 
K ~ 3 GHz yields t c ~ 100 ps, the duration of the single- 
shot readout pulse should be not longer than several tens of 
picosecond. Thus, the weaker interbit coupling, e.g., lowered 
to hundreds of KHz, is required for the current SET tech- 
nique, whose response time is usually hundreds of nanosec- 
ond EH. 

In summary, based on conditionally manipulating the se- 
lected qubits, we have shown how to engineer the macro- 
scopic quantum entanglement of Josephson qubits with fixed- 
couplings. Our proposal allows to deterministic ally pre- 
pare three-qubit GHZ entangled states and allows a macro- 
scopic test of the contradiction between the noncommutativ- 
ity of quantum mechanics and the commutativity of classical 
physics. 
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